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Here, we study the effects of stochastic nuclear motions on the electron transport in doped polymer 
fibers assuming the conducting state of the material. We treat conducting polymers as granular 
metals and apply the quantum theory of conduction in mesoscopic systems to describe the electron 
transport between metalliclike granules. To analyze the effects of nuclear motions we mimic them by 
a phonon bath, and we include electron-phonon interactions in consideration. Our results show that 
the phonon bath plays a crucial part in the intergrain electron transport at moderately low and room 
temperatures, suppressing the original intermediate state for the resonance electron tunneling, and 
producing new states which support the electron transport. Also, the temperature dependence of 
magnitudes of the peaks in the electron transmission corresponding to these new states is analyzed. 

PACS numbers: 72.15.Gd,71.18.+y 



I. INTRODUCTION 

Currently, electron transport in conducting polymer 
nanofibers is a subject of intense interest because such 
nanofibers are expected to be used in building up na- 
noelectronic devices. Doped polyacetylen or polyanilinc- 
polyethylene oxide fibers as well as polypyrrole nanotubes 
could behave as conductors, and their conductivity signif- 
icantly increases upon doping. Charge transport mecha- 
nisms in the conducting polymers were intensively stud- 
ied starting from the very discovery of these materials 

It is known that chemically doped conducting poly- 
mers are very inhomogeneous. In some regions polymer 
chains are disorderly arranged forming an amorphous 
poorly conducting substance. In other places the chains 
are ordered and densely packed [!, Q . The correspond- 
ing regions behave as metalliclike grains embedded in 
the disordered environment. The key point in studies 
of the electron transport in conducting polymers is to 
elucidate the nature and physical mechanisms of the in- 
tergrain transport. Prigodin and Epstein suggest that 
the intergrain transport mostly occurs due to the elec- 
tron tunneling between the grains through intermediate 
resonance states on the polymer chains connecting them. 
This approach was employed to build up the theory of 
electron transport in polyaniline based nanofibers [B[ pro- 
viding a good agreement with the previous transport ex- 
periments Q. Assuming the electron tunneling through 
the intermediate state to be the predominant mechanism 
for the intergrain transport, we see a similarity in elec- 
tron transport mechanisms in conducting polymers and 
those in molecules connecting metal leads. In the case of 
polymers metalliclike domains take on part of the leads, 
and the molecular bridge in between is reduced to a single 
intermediate site. 



An important issue that arises in the analysis of the 
intergrain electron transport is the influence of nuclear 
motions in the medium between the grains. The sim- 
ilar issue was analyzed while developing the theory of 
conduction through macromolecules. It was shown that 
the environment acts as a source of incoherence for the 
tunneling electron 0, H, Ij jjlfjl l . A lso, it can give rise to 
some extra electron states 12[ • In the previous works 
[S [l3|, [3| the effects of stochastic nuclear motions were 
analyzed by introducing the coupling of the intermidiate 
state between the grains to a phonon bath including a 
large amount of harmonic oscillations. It was shown that 
the direct coupling of the bridge to the dissipative reser- 
voir (phonon bath) brings an inelastic component to the 
intergrain transport. This erodes the peak in the elec- 
tron transmission corresponding to the electron tunnel- 
ing through the bridge. Nevertheless, estimations made 
in the recent work [15| give grounds to believe that the 
peak could be still rather well pronounced under typical 
conditions of experiments on the electrical characteriza- 
tion of conducting polymer nanofibers. 

However, the studies of dissipative effects in the in- 
tergrain electron transport in conducting polymers may 
not be restricted with the plain assumption of the direct 
coupling of the bridge site to the phonon bath. Other 
scenarios can occur. In particular, analyzing the electron 
transport in polymers, as well as in molecules, one must 
keep in mind that besides the bridge sites there always 
exist other nearby sites. In some cases the presence of 
such sites may strongly influence the effects of stochastic 
nuclear motions on the characteristics of electron trans- 
port. This may happen when the nearby sites somehow 
"screen" the bridge sites from direct interactions with the 
phonon bath. Here, we elucidate some effects which could 
appear in the electron transport in conducting polymer 
fibers in the case of such indirect coupling of the bridge 
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state to the phonon bath. 



II. MODEL AND RESULTS 

We mimic the effects of the environment by assum- 
ing that the side chain is attached to the bridge, and 
this chain is affected by the phonon bath. This model re- 
sembles those used to analyze electron transport through 
macromolecules [H, E3] • The side chain is introduced 
to the model to screen the resonance state making it more 
stable against the effect of the phonon bath. As in the 
previous papers 13, 12 1 we assume that electrons cannot 
hop along the side chain, so it may be reduced to a single 
site attached to the resonance site (the bridge). 

Within the adopted model the Hamiltonian takes the 
form: 



H — H e i + H p h + H e - p h + Hl + Hr (1) 

where the Hamiltonian H e i describes the bridge with the 
attached side site: 



H el = eb + b + ec + c - w(b + c + h.c), 



(2) 



Hl,r are self-energy terms describing the coupling of 
the left (L) and right (R) grain to the bridge, and e 
and e are the on-site energies of the bridge site and the 
attached side site, respectively. The factor w is the hop- 
ping integral between the bridge an the attached site. 

The remaining terms in the expression (1) describe 
the phonon bath and the electron-phonon interaction, 
namely: 

H ph + H e _ ph = Y / n a B+B a +Y / ^c + c{B a + B+) (3) 

a a 

where Q are the frequencies of phonons belonging to 
the bath and A Q denote the electron-phonon coupling 
strengths. Performing the unitary transformation H = 
e s He~ s with the generator S [la]: 



5 = E^r c+c (^-i^), 



(4) 



the coupling to the bath is eliminated, and we obtain: 

H = H e jf + Hph (5) 

The effective Hamiltonian for the bridge within the 
adopted model takes on the form: 

H eff = H el +H L + H R . (6) 

Here, the Hamiltonian H e i may be presented as follows: 

H e i = Hb + H c + Hb- C (7) 

where again H^ = eb + b corresponds to the bridge site, 
H c is the Hamiltonian of the side chain, and i?&_ c de- 
scribes the interaction between them. The terms H r and 



Hb-c are modified due to the coupling of the side chain 
to the bath. Keeping only one site in the side chain we 
have: 



H c = (e - S)c+c, 

H b -c = -wx(b + c + h.c). 
The quantities 5 and \ are given by [101 ]: 



X = exp 



(8) 



(9) 



(10) 



The lowest order in the Fourier transform for the 
Green's function associated with the Hamiltonian (6) 
reads: 

G-\E) =E-e-E L (E)- E R (E) - w 2 P(E). (11) 

The first four terms in this expression represent the in- 
versed Green's function for the resonance site coupled to 
the two grains where Til,r(E) are complex self-energy 
terms. The last term represents the effect of the environ- 
ment and has the form [Ioj |: 



P(E) 



dtexp[it(E 



S + i0+)} 



x {(1 - /) exp[-F(t)] +/exp[-F(-t)]} (12) 



with exp[— F(t)] being a dynamic bath correlation func- 
tion, and / taking on values 1 and when the attached 
site is occupied and empty, respectively. 

Characterizing the phonon bath with a continuous 
spectral density J(u>) given by [lij |: 



Jo ( — ) exp 



to 



(13) 



one may write out the following expressions for the func- 
tions F{f) and 5 : 



F(t) 



J{u) 



1 



2[1 - cos(c^)] 
exp(w/6') — 1 



S 



-J(lo) = J 



(14) 
(15) 



where 9 is the temperature expressed in the units of 
energy. 

Within the short time scale (oj c t <§; 1) the function 
F(t) could be presented in the form: 



F(t) = — \ - ln[l + (uj c t) 2 } + i arctan(w c i) + Kit) 

UJ r 2 



where 



K(t) = 9 2 t 2 ( [2- 1 — + i)=e 2 t 2 J2 



(16) 
(17) 
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Here, f (2;9/uj c + l) is the Riemann £ function [lg ]. 
The asymptotic expression for K(t) depends on the rela- 
tion between two parameters, namely, the temperature 9 
and the cut-off frequency lj c characterizing the thermal 
relaxation rate of the phonon bath. Assuming 9 uj c 
and estimating the sum of the series in Eq. (12) and we 
get: 



K(t) « —(u c tf 



In the opposite limit lu c ^> 9 we obtain: 



Kit) 



7T 

~6 



{9tf 



(18) 



(19) 



Also, we may roughly estimate K(t) within the inter- 
mediate range. Taking 9 sa uj c we arrive at the approx- 
imation K(t) w a 2 (9t) 2 where a 2 is a dimensionless 
constant of the order of unity. Correspondingly, within 
the short time scale we can omit the first term in the 
expression (16), and we get: 



F(t) » — {iu c t + K{t)} 

OJc 



(20) 



where K(t) is given by either Eq. (18) or Eq. (19) 
depending on the relation between uj c and 9. 

Within the long time scale ui c t ^S> 1, and provided 
that temperatures are not very low {9 3> lu c ), we may 
present the function K (t) as: 



K (t) = 29t 



dz 

-^-(1 — cos z)e 



■z/uj c t 



w9t. 



(21) 



Now, the term K(t) is the greatest addend in the ex- 
pression for F(t), so the latter could be approximated 
as: F{t) w Tt9tJo/uJc- The same approximation holds 
within the low temperature limit when w c >fl>t _1 . 

Using the asymptotic expression (20), we may calcu- 
late the contribution to P(E) coming from the short 
time scale (cu c t <C 1). It has the form: 



Pi(£;)= -W^ exp 





jl + $ 




} 


[ 4J # 




. 2 v / Jo^ _ 





(22) 



where <&(z) is the probability integral. When both ui c 
and 9 have the same order of magnitude the expression 
for P{E) still holds the form (22). At 9 < uj c , the 
temperature 9 in the expression (22) is to be replaced 
by uj c . We remark that under the assumption 9 uj c 
the function P\{E) does not depend on the cut-off fre- 
quency uj Ci whereas at to c 3> 9 it does not depend on 
temperature. The long time (co c t 3> 1) contribution to 
P{E) could be similarly estimated as follows: 



P 2 (E) = 



1 



E - I + Jo 



rJ 9/u 



(23) 




FIG. 1: Energy dependence of the real (top panel) and imag- 
inary (bottom panel) parts of P(E). The curves are plotted 
assuming J = 20rneV, w = WOmeV, I = 0, T = lOOif 
(dash lines), and T — 300K (solid lines) at 6 3> lo c . 



Comparing these expressions (22) and (23) we see that 
the ratio of the peak values of P 2 (E) and Pi (E) is of 
the order of (uj 2 /J 9) 1 ^ 2 . Therefore the term P\{E) pre- 
dominates over P2(E) when the temperatures are mod- 
erately high (lo c < 9) and the electron-phonon interac- 
tion is not too weak Jq/uj c ~ 1. Usually, experiments 
on the electrical characterization of conducting polymer 
nanofibers are carried out at T ~ 100 -j- 300K, so in fur- 
ther analysis we assume that iuj 2 / J^ff) 1 ^ 2 -C 1, and we 
omit the term P2(E) from the following consideration 
using the approximation P(E) w Pi(E). As shown in 
the Fig. 1, the imaginary part of P(E) exhibits a dip 
around E — e and the width of the latter is determined 
by the product of the temperature 9 (or lo c ) and the 
constant Jo characterizing the strength of the electron- 
phonon interaction. When either factor increases, the dip 
becomes broader and its magnitude reduces. 

The presence of P(E) gives rise to very significant 
changes in the behavior of the Green's function given by 
the Eq. (11). Using the flat band approximation for the 
self-energy corrections namely: Sl,_r. = — i^-L.R 

where A.l,r are constants of the dimensions of energy, 
and neglecting for a moment all imaginary terms in the 
Eq. (11), we find that two extra poles of the Green's func- 
tion emerge due to the environment. When the phonon 
bath is detached (Jo = 0) the separation between the 
poles is determined with the hopping integral w. Due to 
the electron coupling to the phonons the positions of the 
poles become shifted. Assuming e = e = and 9 uj c 
we get: 



E = ±2^/j 9\\n(2J 9/w 2 )\. 



(24) 
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The poles correspond to extra electron states which oc- 
cur due to electrons coupling to the environment. These 
new states are revealed in the structure of the electron 
transmission T(E) which reads [2(j |: 



T(E) = 4A L A R \G(E)f 



(25) 



The structure of T(E) is shown in the Fig. 2. Two 
peaks in the transmission are associated with the envi- 
ronment induced electronic states. Their positions and 
heights depend on the temperature and on the coupling 
strengths Jo, and w. The important feature in the elec- 
tron transmission is the absence of the peak associated 
with the resonance state between the grains (the bridge 
site) itself. This happens due to the strong suppression of 
the latter by the effects of the environment. Technically, 
this peak is damped for it is located at E = where the 
imaginary part of P(E) reachs its maximum in magni- 
tude. To provide the damping of the original resonance 
the contribution from the environment (including the side 
chain attached to the bridge) to the Green's function (11) 
must exceed the terms £l,.r describing the effect of the 
grains. This occurs when the inequality 



A < w 2 /, 



(26) 



is satisfied. When the coupling of the bridge to the at- 
tached side site is weak, the influence of the environment 
slackens and the original peak associated with the bridge 
at E = e may emerge. At the same time the features 
originating from the environment induced status become 
small compared to this peak. 

So, the effects of the environment may lead to the 
damping of the original resonance state for the electron 
tunneling between the metallic islands in the polymer 
fiber. Instead, two environment induced states appear 
to serve as intermediate states for the electron trans- 
port. Similar effects were recently investigated in the 
electron transport through DNA macromolecules fioj ] . 
and it was shown that low biased current-voltage charac- 
teristics may be noticeably changed due to the occurence 
of the phonon bath induced electron states. Due to some 
particular features of conducting polymers, these effects 
may be revealed in polymer nanofiber current-voltage 
curves at significantly higher bias voltage, as discussed 
below. 

Realistic polymer nanofibers have diameters within the 
range 204- 100 nm, and lengths of the order of a few mi- 
crons. This is much greater than the typical size of both 
metalliclike grains and intergrain separations which take 
on values ~ 5 4 10 nm 21|. Therefore, we may treat a 
nanofiber as a set of parallel working channels, any sin- 
gle channel being a sequence of grains connected with the 
resonance polymeric chains. The net current in the fiber 
is the sum of currents flowing in these channels, and the 
voltage V applied across the whole fiber is distributed 
among sequental pairs of grains along a single channel. 
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FIG. 2: The renormalized electron transmission function 
vs energy. The curves are plotted for e = e = 0, w = 
WOmeV, T = WOK (dash lines) and T = 300K (solid lines). 
The constant Jo equals 20meV (top panel) and 50meV 
(bottom panel). 



So, the voltage AV applied across two adjacent grains 
could be estimated as AV^ ~ VL/Lq where L is the av- 
erage separation between the grains, and Lq is the fiber 
length. The ratio AV jV may take on values of the order 
of 10 _1 -^10 -2 . Experiments on the electrical character- 
ization of the polymer fibers are usually carried out at 
moderately high temperatures (T 4 300K) , so it seems 
likely that 9 > u> c . Assuming that w ~ lOOme^, and 
Jo ~ 20450meV we estimate the separation between the 
environment induced peaks in the electron transmission 
as 120 4l70mey (see Fig. 2). This estimate is close to 
eAV when V takes on values up to 2 4 3 volts. So, the 
environment induced peaks in the electron transmission 
determine the shape of the current-voltage curves even at 
reasonably high values of the bias voltage applied across 
the fiber. 

In calculations of the current we employ the formula 



1=^1 dET(E)[f 1 (E) - f 2 (E)}. 



(27) 



Here, n is the number of the working channels in the 
fiber, f 12(E) are Fermi functions taken with the differ- 
ent contact chemical potentials [ii^. for the grains. The 
chemical potentials differ due to the applied bias voltage 
A^ : 

iii =E F + (1- rj)eAV; [x 2 = E F - rjeAV. (28) 

The parameter 77 characterizes how the voltage AV^ is 
divided between the grains; Ep is the equilibrium Fermi 
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FIG. 3: The current-voltage characteristics (nA-V) plotted 
for n=10, Ac = A R = 0.5meV, rj = 1/2, w = WOmeV, 
Jo = 20meV (left panel), and Jo = 50meV (right panel) at 
T = WOK (dash lines) and T = 300^ (solid lines). 



energy of the system including the pair of grains and 
the resonance site in between, and T(E) is the electron 
transmission function given by Eq. (25). 

The resulting current- voltage characteristics are shown 
in the Fig. 3. The current takes on low values (~ InA) 
because the coupling of the grains to the intermediate 
state is weak due to comparatively large distances be- 
tween the grains 0] ■ Consequently, A.l,r take on values 
much smaller than those typical for electron transport 
through molecules placed between metal leads. The I—V 
curves exhibit a nonlinear shape even at room tempera- 
ture despite the fact that the original state for the res- 
onance tunneling is completely suppressed. This occurs 
for the intergrain transport is supported by new phonon 
induced electron states. In this work we did not take a 
task of providing a quantitative agreement between the 
theory and experimental results obtained on conducting 
polymer fibers. However, in general, our calculated I—V 
curves agree with the reported experiments Q- 

It is worthwhile to discuss the temperature dependence 
of the peak value of the electron transmission which 
follows from the present results. Using the expression 
(11) for the Green's function and the expression (22) for 
P{E) we may conclude that at low temperatures the 
transmission accepts small values, and exhibits rather 
weak temperature dependence. At higher temperatures 
(T ~ 100JQ the transmission increases fast as the tem- 
perature rises and then it reduces as the temperature fur- 
ther increases. The peak in the transmission is associated 
with the most favorable conditions for the environment 
induced states to exist when all remaining parameters 
(such as Jo and w) are fixed. At high temperatures the 
peaks associated with the environment induced states are 
washed out, as usual. This is shown in the Fig. 4. Also, 
basing on the results of the previous works @, [IH we 
may analyze the temperature dependence of the electron 
transmission function assuming that the bridge between 




FIG. 4: Arrhenius plot of the peak value of the electron trans- 
mission function for Jo = 20meV, w = WOmeV, Al = Ar = 
0.5meV; 6 = 8.6meV (T = 100JC). Solid line is plotted as- 
suming the indirect coupling of the bridge to the phonon bath 
(Eqs. (11), (25)). Dashed line is plotted assuming that the 
bridge is directly coupled to the phonons (Eq. (29)). 



two adjacent grains is directly coupled to the phonon 
bath. In this case the transmission peak value may be 
presented in the form [l5| : 



where 



T(E) = 1 - 



16J 



2(1 + v/TiV) 2 



:C 2 



l 



(29) 



(30) 



LU C (A L +A R ) 

The temperature dependence of the transmission given 
by the Eq. (29) is also shown in the Fig. 4. Both curves 
are plotted at the same value of the electron-phonon cou- 
pling strength Jo . Comparing them we conclude that at 
higher temperatures the curves appear to differ. While 
the temperature rises, we observe a peak in the elec- 
tron transmission assuming the indirect coupling of the 
bridge to the phonon bath, and we see the transmission 
to monotonically decrease when we consider the bridge 
directly coupled to the bath. Correspondingly, we may 
expect qualitative diversities to appear in the tempera- 
ture dependencies of the current, as well. These diver- 
sities originate from the difference in the effects of envi- 
ronment on the intergrain electron transport in the cases 
of direct and indirect coupling of the bridge site to the 
phonon reservoir. When the bridge is directly coupled to 
the bath the stochastic motions in the environment only 
cause the peak in the electron transmission to be eroded, 
and the higher is the temperature the less distinguish- 
able is the peak. However, when the bridge is screened 
from the direct coupling with the phonons due to the 
presence of the nearby sites, the stochastic nuclear mo- 
tions in the medium between the grains (especially those 
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in the resonance chain) may take a very different part in 
the electron transport in conducting polymers at moder- 
ately low and room temperatures. Due to their influence 
the original intermediate state for the resonance tunnel- 
ing may be completely suppresed but new environments 
induced states may appear providing the electron trans- 
port through the polymer nanofibers. 

III. CONCLUSION 

Finally, studies of the electron transport in conduct- 
ing polymers are not completed so far. There are solid 
grounds to expect significant dissipative effects in the in- 
tergrain transport at moderately high temperatures. As 
shown in the previous studies of the electron transport 
through molecules, various dissipative effects may occur 
depending on characteristic features in the interaction 
of a propagating electron with the environment. Among 
these features we may single out the character of the elec- 
tron coupling to the dissipative reservoir (phonon bath) 
as a very significant factor. It is likely that in realistic 
conducting polymers both direct and indirect coupling 
of the intermediate state (the bridge) to the environment 
may occur. We need more experimental data concern- 
ing transport in these materials to determine which kind 
of coupling would prevail in a particular material under 
particular experimental conditions. The model adopted 
in the present work predicts that the temperature de- 
pendence of the electron transmission through the envi- 
ronment induced states crucially differs from that which 
is typical for the electron transport through the original 
bridge state directly coupled to the dephasing reservoir. 
So, the present results give means to experimentaly study 
distinctive features of the dissipative effects in the elec- 
tron transport through conducting polymer nanofibers. 
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